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Faster quantum searching with almost arbitrary operators 
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Grovers search algorithm drives a quantum system from an initial state |s) to a desired final state 
|t) by using selective phase inversions of these two states. In 01, we studied a generalization of 
Grovers algorithm which relaxes the assumption of the efficient implementation of A, the selective 
phase inversion of the initial state, also known as diffusion operator. This assumption is known 
to become a serious handicap in cases of physical interest @Ta|. Our general search algorithm 
works with almost arbitrary diffusion operator Ds with only restriction of having |s) as one of 
its eigenstates. The price that we pay for using arbitrary operator is an increase in the number 
of oracle queries by a factor of 0{B), where i? is a characteristic of the eigenspectrum of Ds 
and it can be large in some situations. Here we show that by using quantum fourier transform, 
we can regain the optimal query complexity of Grovers algorithm without losing the freedom of 
using arbitrary diffusion operators for quantum searching. However, the total number of operators 
required by algorithm is still 0{B) times more than that of Grovers algorithm. So our algorithm 
offers advantage only if oracle operator is computationally more expensive than diffusion operator, 
which is true in most search problems. 

PACS numbers: 03.67.Ac 


I. INTRODUCTION 

Suppose we have a quantum computer initially in a 
source state |s) and we want to evolve it to a final state 
|t), say the target state. One way to do this is famous 
Grovers algorithm which assumes that we can effi¬ 
ciently implement the selective phase inversion operators, 
Is and It, of these two states. Grovers algorithm keeps 
on iterating the Grovers search operator Q = Iglt on the 
source state |s) to get the target state |t). The required 
number of iterations is Tr/dcr, where a = |(t|s)|. 

While using Grovers algorithm for the search problem, 
we choose |s) to be the uniform superposition of all N ba¬ 
sis states to be searched i.e. |s) = \i)/^fN. In case of 

a unique solution to the search problem, the target state 
|t) is a unique basis state and a = |(t|s)| = I/VN. Thus 
Grovers algorithm outputs a solution in just 0{y/N) time 
steps which is quadratically faster than classical search 
algorithms taking 0{N) time steps. 

Grovers algorithm is proved to be strictly optimal [ 1 ] 
and the assumption of efficient implementation of Ig and 
It is justified in many situations. The oracle operator It 
always needs an oracle query for implementation whereas 
implementation of the diffusion operator Is is dictated by 
physical constraints and sometimes it may become a seri¬ 
ous problem. For example, in the case of two-dimensional 
spatial search 0 , each implementation of Ig takes 
time steps and hence the total time complexity of Grovers 
algorithm becomes \/]V x \/]V = N time steps, which is 
no better than classical algorithms. 

Earlier, several attempts have been done to study 
quantum search algorithms with more general diffusion 
operators in place of Is- For example, Kato studied the 
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case when Ig is replaced by an operator made up of only 
single qubit gates Q- Ambainis studied the case when 
Is is rralaced W a real operator with |s) as its eigen¬ 
state [ 3 . In [l|, we presented a unified framework for 
all such generalizations and we studied the case when Ig 
is replaced by any arbitrary operator Dg with the only 
restriction of having |s) as one of its eigenstates. This re¬ 
striction seems to be more or less justified as the search 
operator should have some special connection with the 
source state. 

To be more concrete, let the normalized eigenspec¬ 
trum of Dg be given by Dg\i) = with \t) as the 

eigenstates and {9t) as the corresponding eigenval¬ 
ues (eigenphases). Since a global phase is irrelevant in 
quantum dynamics, we choose Ds\s) = |s), i.e. Oi^g = 0. 
In 0, we studied the iteration of general search opera¬ 
tor S = Dgit on |s) by analyzing the eigenspectrum of 
S. We found that the performance of quantum search 
algorithm depends upon the quantities Ai and A 2 , where 

At, = j2\m\^oot^^ (1) 

i^s 

is the moment of cot ^ with respect to the distribu¬ 
tion |(^|t)P over all if^s. 

We found that Al should be very close to zero for a 
successful quantum search. We also presented an algo¬ 
rithm in Section IV.A of 0 which uses an ancilla qubit 
to control the applications of Dg and D\ in a clever way 
so that effectively Ai becomes zero. In this paper, we 
restrict ourselves to the case Ai = 0. In case, it is not 
so, we can always use just-mentioned algorithm to make 
it so. If Al = 0, we get the target state using {'K/Aa)B^ 
oracle queries where B = \/\ + A 2 . In Section IV.B of 
0 , it was shown that by controlling the operators us¬ 
ing an ancilla qubit, we can get a faster algorithm which 
uses only ( 3 -\/ 3 / 2 )( 7 r/ 4 a)i 3 oracle queries. But even now. 
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this algorithm is slow compared to the optimal Grovers 
algorithm by a factor of 0{B). 

In this paper, essentially, we show that by using quan¬ 
tum fourier transform, we can manipulate the opera¬ 
tion of Dg cleverly to design a new operator for which 
B = 0(1). Thus we can achieve the optimal perfor¬ 
mance of Grovers algorithm with almost any operator 
Dg. In next section, we present a brief review of the re¬ 
sults of our analysis of general quantum search algorithm 
presented in [I|. In Section III, we present our new al¬ 
gorithm with optimal performance. Then we conclude in 
Section IV. 


II. GENERAL QUANTUM SEARCH: A BRIEF 
REVIEW 


Here we briefly discuss the dynamics of general quan¬ 
tum search algorithm presented in [T|. We skip the de¬ 
tails for which readers are referred to the original paper. 
This algorithm iterates the operator S — Dgif on |s) to 
take it close to |t). Here Dg is as defined earlier and if 
is the selective phase rotation of target state by angle of 
(j). Without any loss of generality, we take ^ = tt so that 
if is the selective phase inversion It of the target state. 
Also, we assume |s) to be a non-degenerate eigenstate 
for simplicity. Let the normalized eigenspectrum of Dg 
be given by Ds\£) = By convention, 9g^g = 0. 

Let other eigenvalues satisfy 

I > 6*111111 > 0, Oie[-Tr,n] (2) 

To study the iteration of 5 on |s), we need to find its 
eigenspectrum. We found that only two eigenstates |A±) 
with the corresponding eigenvalues of S are relevant 
for our algorithm. Under the assumption, |A±| <C 0min, 
the initial state |s) is almost completely spanned by two 
eigenstates |A±). The eigenvalues A± are given by 

A± = ±^(tan?7)±^ cot 277 =^^. (3) 

where 

H = Ap = ^|(£|t)pcot^>| . (4) 

i^s 


As a = \{s\t)\ <C 1, J2e^g l(^|6)P is very close to 1 and 
hence above equation gives 



( 5 ) 

We only consider the case when Ai = 0 as only then 
a successful quantum search is possible. If Ai ^ 0, we 
can always use the algorithm presented in Section IV.A 
to design a new operator using Dg and Dl for which 
Ai = 0. In this case, |3] indicates that 


TT 2a 

^= 4 ’ = 


With T] = n/I and (j) = tt, Eq. (23) and (24) of [l| 
gives us the initial state |s) and the effect of iterating S 
on |s) in terms of two relevant eigenstates. We have 

|s) = -*/y2[e*^+/2|A+)-e*^-/2|A_)], (7) 

and 

=-*/V2[e*«'^+|A+)-e"«'"-|A_)], (8) 

where q' = q + 

For q = q^ Ki 7r/2|A±| = irB/Aa, the state is 

very close to the state |?«) given by 

S‘>-\s) = \w) = l/V2i\X+) + \X.)). (9) 

As shown in [ij, we have |(t|rc)| = 1/5. So we get the 
target state with a probability of 1/5^ after 7r5/4a it¬ 
erations of S on the initial state |s). As each application 
of S needs one query to implement It , we need a total of 
(7r/4a)5'^ oracle queries to get the target state. Section 
IV.B of [l[ presents an algorithm which improves the 
query complexity to 0{'!tB jAa) by controlling the appli¬ 
cations of Dg using an ancilla qubit. 

However, we still need 0{B) times more queries com¬ 
pared to the optimal Grovers algorithm. We can see from 
|4]that A 2 < l/6*min so for dmin 'C 1, the only upper bound 
that we have on B is 1/dmin, which can in general be ar¬ 
bitrarily large. Thus we need some trick to reduce the 
value of B. One way is to use quantum fourier transform 
to reliably distinguish the eigenstates of Dg using Phase 
Estimation algorithm. As 9i=s = 0 and \9i^s\ > 0minj we 
need O(l/0min) applications of Dg to achieve such a reli¬ 
able distinction after which, we can selectively invert the 
phase of |s) state to implement Ig. Such a scheme will 
take O(l/a0niin) applications of Dg and 0(l/a) oracle 
queries. 

We point out that 1/dmin is just an upper bound for 
B hence above-mentioned simple scheme may not be 
more efficient than general quantum search algorithm 
with time complexity 0{B/a). For example, in case 
of two-dimensional spatial search, l/0min = 0{'/N) but 
B — C)(-\/ln N) so B is much smaller than l/dmin- In next 
section, we show that the quantum fourier transform can 
be used in a more clever way to get a successful quan¬ 
tum searching using only 0{B/a) applications of Dg and 
0(1 /q!) oracle queries. 


III. FASTER ALGORITHM 


The new algorithm basically reduces the effective value 
of 5 to 0(1). To get the basic idea, first consider a 
simple case when the operator Dg is iterated r times to 
get a new operator with its eigenspectrum given by 
D^g\t) = As given byO the corresponding value 

of {BrY is 


{Brf = V 

l^S 


\m\^ 


Ai = 0 


( 6 ) 


2 


( 10 ) 
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As long as 0^/2 <C l/r, we have sin(r0£/2) Ri rsin(0£/2) 
and so ~ Bjr. Thus can be made 0(1) by choos¬ 
ing r = 0{B). But this simple scheme works only if 
9((2 1/r for all i, not in the general case. In fact, we 

can see that if 9^ is sufficiently close to 2mTlr for any 
integer n, then Br will diverge. 

We now show that the quantum fourier transform can 
be used in more general case to make B = 0(1). Ba¬ 
sically we achieve some level of distinction among the 
eigenstates of Dg and then to effectively nullify the con¬ 
tribution of those eigenstates for which 9^/2 is not much 
smaller than 1/r. 

Let Hn denote the Hilbert space of our main quantum 
system which we want to evolve from |s) to \t) state. 
We attach an ancilla quantum system of m qubits to 
our main system and let Hm {M = 2 "*) denote the 
corresponding Hilbert space with its basis states |j), 
j G {0,1,..., 2"* — 1}. This ancilla system will be used 
for quantum fourier transform as done in Phase Esti¬ 
mation algorithm. We work in the joint Hilbert space 
Ti = T~Lm (SiBn- 

A. Phase Estimation Algorithm 


As the main quantum system is in an eigenstate |£) of 
Dg, we get 

= i E (15) 

j=0 

The third operator P 3 = IF 'Si In applies a quan¬ 
tum fourier transform on the ancilla quantum system but 
leaves the main system unchanged. The action of T in 
a 2™-dimensional Hilbert space on each basis state \j) is 
given by 

2--1 

E exp{2Tnkj/2"^) \j). (16) 

k=0 

where k is the numerical value of the binary number rep¬ 
resented by the bit string encoded by \k) state. For ex¬ 
ample, for \k) = | 0 ), all qubits are in | 0 ) state and hence 
k = 0. 

As suggested by Uni before applying T, the ancilla 
quantum system is in the state (1/2™/^) J2j So, 

using m we get 


We first consider the operator V corresponding to the 
phase estimation algorithm. Let the initial state of our 
main quantum system be |£), an eigenstate of Dg, and 
let the initial state of our ancilla quantum system be | 0 ) 
in which all m qubits are in |0) state. So the initial state 
for PEA is 

| 0 ,£) = | 6 )m®K)^, ( 11 ) 

where the subscripts M and N denote the corresponding 
Hilbert space of quantum states. We omit these sub¬ 
scripts for simplicity. 

The phase estimation algorithm operator P is a suc¬ 
cessive application of three operators on the state | 0 ,^), 

i.e. 

V = {TS1n){cjDI){WsIn)\0J). (12) 

The first operator Pi = W S In applies Walsh- 
Hadamard transform on the ancilla system and leaves 
the main system unchanged. Thus the ancilla system is 
transformed to an uniform superposition of all m basis 
states and we get 


p\6j) = mi). ( 17 ) 


where 


^ E ea;p[i(27rjfc/2”"j6»£)]|fc). (18) 

k,j=0 

This is the standard output state of Phase Estimation 
algorithm and has been analyz ed q uite well in literature 
(for example, see Sec. 5.2.1 of [l3)- We note that 

2 "*-! 

^ E {exp[t{27Tk/2”^ - 9i)]y, (19) 

j=o 


which is the sum of a geometric series and after little 
calculation, we get 


\{km = 


1 sin[7rfc2'" ^ 9 (\ 
^sin[(7rfc2""-i6»<!)/2™] ‘ 


( 20 ) 


\i)v). ( 13 ) 

j=0 

The second operator P 2 = cjDi is a controlled appli¬ 
cation of D/g operator on the main quantum system, i.e. 
it applies j iterations of Dg on the main quantum system 
if and only if the ancilla quantum system is in the \j) 
state. 

^27’i|6,£) = ^ (14) 

1=0 


For the purpose of our algorithm, what matters is the 
overlap of \9(,) state with \k = 0) state for which k = 0 
and above equation gives us 

_ 1 sin[ 2 — 10 ,] 

2^ sin[ 6 »f/ 2 ] ■ ^ ’ 

For i = s, 9i = 0 and in the limit 9i —>■ 0, the R.H.S. of 
above equation becomes 1 so we have |0s) = |0). So if the 
main quantum system is in |s) state, then the operator 
P leaves the entire quantum system unchanged. 


I( 6 | 0 i)| 
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B. New operator 


So the corresponding contribution to (-B)^ is 


The new operator that we design using Dg as basic 
operator is given hy V = VCV\ where C is a diagonal 
operator given by 


C = (c,^6(-l^))(co^r)- (22) 

First, C applies the operator which performs r iter¬ 
ations of the operator Dg on the main quantum system, 
if and only if the ancilla quantum system is in jO) state. 
ThenC applies the operator c^^q(—Iat) which inverts the 
phase of all basis states \j) of the ancilla quantum sys¬ 
tem except the \j = 0) state. It is easy to check that 
the eigenstates and the corresponding eigenvalues of the 
operator C are given by 

|j = 6)|^) ; exp[ir0i] 

\j^0)\i) ■ exp[nr] . (23) 

As V — VCP^, we find that the eigenstates \t} and the 
corresponding eigenvalues of our new operator D are 
given by 

^(|j=0)K)) ; explirOi] (24) 

^(|j7^0)K)) ; exp[nT] . 


We choose 

r = 2"* . (26) 

Then it is easy to check that the numerator of one term in 
[25] exactly cancels the denominator of another term. It is 
this cancellation which effectively nullifies the contribu¬ 
tion of those eigenstates which otherwise had a potential 
to make B diverge as we discussed earlier. And as we 
have shown that such a cancellation can be achieved eas¬ 
ily through quantum fourier transform. After this can¬ 
cellation, we get 

1 V- \{W ^ 

^2 — 22m ^2 gj 2 Oe ~ 22m ' (2'^) 

^ i^s T ^ 

And so, we get 

(5)2 = Sj -I- E 2 < 1 + ■ (28) 


The eigenstate corresponding to = 0 is the effective 
source state |s) and above equation gives |s) = B(|0)|s)). 
As discussed in the previous subsection, when the main 
quantum system is in |s) state, then the operator V leaves 
the entire quantum system unchanged. So we have |s) = 
|6)|s). 

Now consider the oracle operator used by our algo¬ 
rithm. In place of simply applying the selective phase 
inversion It on the main quantum system, our new al¬ 
gorithm applies It on the main system if and only if the 
ancilla system is in |0) state. By doing this, basically we 
apply Ig j, the selective phase inversion of |0)|t) state and 
our effective target state is \t) = |0)|t). 

With all these details of new operators, the value of 
B for T> can be found. We use the expression jS] i.e. 
(5)^ = 'T^i^g{W\'t)\/ sin y)^- We separate this sum into 
two parts: Si and S 2 where Si is due to the eigenstates 
V{\j ^ 0)1^)) having eigenvalue e™, i.e. Oi = tt. As 
sin^ ^ = 1 for all such eigenstates and the sum I W P 
is always less than 1 for any set of \tj, so we see that 
Si < 1. 

To find S 2 , we note that for remaining eigenstates, the 
eigenvalues are and the values of |(I'|i)| are given by 

|(0|(t|B(|6)|£)) 


which, using is equal to 


\m\moi) 


1 sinp^-^e^j 

2™ sin [0^/2] 


\m\. 


Thus, by choosing m = log 2 B, we can make B = \/2 = 
0(1) and we can get the target state using only 0(7r/4Q;) 
oracle queries, which is within a constant factor of the 
optimal performance by Grovers algorithm. 

However, we note that implementation of the new op¬ 
erator T> requires 3B applications of the basic operator 
Dg, one each for V, 'P\ and C. And as we need 0(7r/4Q:) 
iterations of the search operator, the total applications 
of Dg required by algorithm is still 0{B/a) as required 
by the original algorithm presented in [l| . The improve¬ 
ment mainly comes because of the reduction in number 
of oracle queries which involves major computational re¬ 
sources in typical search problems. 


C. Algorithm 

Above analysis suggests the following algorithm 

(1) We attach an ancilla quantum system of m = log 2 B 
qubits to our main quantum system. 

(2) We choose the initial state of all qubits of our ancilla 
system to be in |0) state. We choose the initial state of 
our main quantum system to be the source state |s). 

(3) We perform 0(1 /q:) iterations of the operator B/g ^ 
on the initial state. Note that Jg ^ is a controlled applica¬ 
tion of It and needs one oracle query for implementation. 
Also, T> needs 0{B) applications of Dg for implementa¬ 
tion. 

(4) We measure the main system, which will be in the 
target state with probability close to 1. 
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IV. DISCUSSION AND CONCLUSION 

We have presented an algorithm to achieve the optimal 
performance of Grovers algorithm in the case of more 
general diffusion operators Dg used in place oi Ig- A 
simple scheme of successive iterations of Dg may not help 
in general cases as some bad eigenstates of Dg may cause 
the parameter B to diverge. What we have basically 
shown is that Quantum Fourier Transform allows us to 
nullify the effect of these bad eigenstates. 

In the example of two-dimensional spatial search a = 
l/\/]V and B = 0(\/lnN). Hence our algorithm can get 
the target state using only 0{'/N) oracle queries in con¬ 
trast to the O(VAlniV) performance of earlier known 


algorithms, see for example 0. However, our algorithm 
will need 0{VN In N) applications of the local operator 
and as each local operator needs one time step f or imple - 
mentation, the total time complexity is still 0{VN In N) 
time steps. But in the cases when oracle query becomes 
more expensive than implementing local operators, our 
algorithm offers an advantage. 

Our algorithm offers a general framework to achieve 
the optimal oracle query performance in general cases. It 
shows that by using general diffusion operators, we dont 
need to compromise with the query complexity of quan¬ 
tum search algorithms. We believe that our algorithm 
can find important applications in search problems. 
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